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Abstract 

We search for decays of a B meson into a neutral D meson and a kaon, with the D me- 
son decaying into K + ir~ir°. This final state can be reached through the b — > c transition 
B~ — ► D°K~ followed by the doubly Cabibbo-suppressed D° — > K + ir~ir°, or the b — > u 
transition B~ —> D°K~ followed by the Cabibbo- favored D° — > K + ir~ir°. The interference 
of these two amplitudes is sensitive to the angle 7 of the unitarity triangle. We present pre- 
liminary results based on 226xl0 6 e+e~ — > T(AS) — > -Bi? events collected with the BABAR 
detector at SLAC. We find no significant evidence for these decays and we set a limit Rads = 
r([K+^-^]pK+)+r([g-7r+7r ( ']jg-) < 0-039 at 95% confidence level, which we translate with a 
Bayesian approach into r B = \A(B~ -> D°K~)/A(B- -> D if")| < 0.185 at 95% confidence 
level. 
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1 Introduction 



Following the discovery of CP violation in B meson decays and the measurement of the angle (3 of 
the unitarity triangle [T] associated with the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix, focus has turned toward the measurements of the other angles a and 7. Following Ref. 2 , 
several methods have been proposed to measure the relative weak phase between the B~ — > D°K~ 
amplitude, proportional to the CKM matrix element V u b (Fig.^), and the B~ —* D°K~ amplitude, 
proportional to V c b- This weak phase, which by definition (7 = arg(— V* b V ud /V* b V cd )) is 7, can be 
measured from the interference that occurs when the D° and the D° decay to common final states. 



U 




Figure 1: Feynman diagrams for the CKM-favored B~ -> D°K~ and the CKM- and color- 
suppressed B~ — > D°K~ decays. 



As an extension of the method proposed in Ref. [Hj, we search for B~ — » [K + tt^tt°]dK~ [3], 
where the CKM-favored B~ —* D°K~ decay, followed by the doubly Cabibbo-suppressed D° — > 
i^ + 7r _ 7r° decay, interferes with the CKM-suppressed _B~ — > D°K~ decay, followed by the Cabibbo- 
favored D — ► i^ + 7r _ 7r° decay. 

In order to reduce the systematic uncertainties, we measure ratios of branching fractions of the 
decay modes of interest in which the two kaons have opposite charge, referred to as "wrong sign" 
events, to the corresponding ones in favored decays, where the kaons have the same charge, referred 
to as "right sign" events. The two ratios we consider, to separate the sensitivity to the suppressed 
rate and the CP violation, are: 

T([K + 7r-TT ] D K-) +T([K-ir+ir°} D K+) 
ads - r ([^+ 7r - 7r 0] D ^+) +r([i^-7r+7r°] D ^-) U 



Aads 



r B + r D + 2r B r£)Ccos7 
r([JC+7r-7r°] D K-) - r([K-TT + TT°] D K + ) 
T([K+tt-tt°} d K-) +T([K-ir+ir°} D K+) 
2rBroS sin 7 /Rads 



where D-mixing effects are neglected, r# = 
parameters take into account the fact that t 



A{B-^D a K~) 



r 2 r, = 1 l^n^^-^-i-^iu : and the C and S 



D — B(D°^K-w+tt°) ' 



A(B-^D U K~) 

re strong phases of the D decays are a function of the 
decay kinematics. Indicating with rrt a point in the Dalitz plane [m 2 KlT , m ^- vr o] ) with [Ad(jH), 5(rft)] 
(\Ad(^), fiffi)]) the absolute value and the strong phase of the D (D) decay amplitudes, and with 
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5 B the strong phase difference between the two interfering B decays, we have 



c _ J Ap{rn)A D (rn) cos(5(rfi) - 5(rft) + 5 B {rft))drn 
yfj \A D (rfi)\ 2 drti ■ J \A D (rft)\ 2 dm 
J A D {rrt)AD(rri) sm(5(Wt) - 5(rft) + 5 B (Wt))drrt 
\A D (rft)\ 2 drft ■ J \A D (nt)\ 2 drn 

Determining the angle 7 from the measurements of Rads an d Aads requires extracting the 
strong phases, for which the available statistics are insufficient. However, the value of r B determines, 
in part, the level of interference between the diagrams of Fig. ^ In most techniques for measuring 
7, high values of r B lead to larger interference and better sensitivity to 7. Thus, r B is a key quantity 
for the extraction of 7 from other measurements in B — > DK decays jH]. In this paper we therefore 
only measure Rads and we constrain r B by exploiting the fact that in Eq. EJ|C| < 1- 

Both the Belle and BABAR collaborations have published similar measurements but in a different 
decay chain (B — > DK with D — > Kir) jSJ. Unlike those measurements, we can take advantage 
of the smaller value of r D , which is r 2 D = (0.214 ± 0.008 ± 0.008)% in D K-ktt decays as 
opposed to rf, = (0.362 ± 0.020 ± 0.027)% in D -> Kir decays. This implies that for a given 
error on Rads, the sensitivity to r B is larger. 



2 Event Reconstruction and Selection 

The results presented in this paper are based on an 226xl0 6 T(4S) BB decays collected between 
1999 and 2004 with the BABAR detector at the PEP-II B Factory at SLAC 9 . In addition, 15.8 ftr 1 
of off-resonance data, with center-of-mass (CM) energy 40 MeV below the T(4S) resonance, is used 
to study backgrounds from continuum events, e + e~ — > qq (q = u,d,s, or c). The BABAR detector 
is described elsewhere jlUj . As far as this study is concerned, charged-particle tracking is provided 
by a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). In addition to 
providing precise spatial hits for tracking, the SVT and DCH also measure the specific ionization 
(dE/dx), which is used for particle identification of low-momentum charged particles. At higher 
momenta (p > 0.7 GeV/c) pions and kaons are identified by Cherenkov radiation detected in a 
ring- imaging device (DIRC). The typical separation between pions and kaons varies from 8a at 
2 GeV/c to 2.5a at 4 GeV/c The position and energy of neutral clusters (photons) are measured 
with an electromagnetic calorimeter (EMC) consisting of 6580 thallium-doped Csl crystals. These 
systems are mounted inside a 1.5-T solenoidal super-conducting magnet. 

The event selection was developed from studies of BB and continuum events simulated with 
Monte Carlo techniques (MC), and of off-resonance data. A large on-resonance data sample of 
B~ — ► D°tt~, D° — ► K~tt + tt° events was used to validate several aspects of the simulation and 
analysis procedure. We refer to this mode as B — > Dn. 

In the reconstruction, both kaon candidates are required to satisfy kaon identification criteria, 
which are based on specific ionization loss measured in the tracking devices and on Cherenkov 
angles in the DIRC and are typically 85% efficient, depending on momentum and polar angle. 
Misidentification rates are at the two percent level. The tt° candidates are reconstructed as pairs 
of photon candidates in the EMC, each with energy larger than 70 MeV and lateral shower profile 
consistent with an electromagnetic deposit, with a total energy greater than 200 MeV, and with 
118.25 < m~ r/ < 145.05 MeV/c 2 . To account for the correlation between the tails in the distribution 
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of the Kirir invariant mass and the tt° candidate mass, we require the difference between the 
two measured masses to be within 32.5 MeV/c 2 of the expected value of 1729.5 MeV/c 2 11 . The 
remaining background from other B^ — ► [/ii^ 71 " ]^^ E! modes is reduced by removing events 
where any /ii/i27r° candidate, with any particle-type assignment except for the signal hypothesis 
for the h\hi pair, is consistent with a D° meson decay. 

After these requirements, the background is mostly due to e + e~ — ► cc events, with c — ► D° — ► 
K + tt^tt° and c — > D — > i'T - . In order to discriminate against them we use a neural network 
(NN) with six quantities that distinguish continuum and BB events: (i) L$ = J2iPi an d (h) 
L2 = cos2 both calculated in the CM frame. Here, pi is the momentum and Oi is the angle 
with respect to the thrust axis of the B candidate of tracks and clusters not used to reconstruct the 
B. (hi) The angle in the CM frame between the thrust axes of the B and of the detected remainder 
of the event, (iv) The polar angle of the B candidate in the CM frame, (v) The distance of closest 
approach between the bachelor track and the trajectory of the D meson. This is consistent with 
zero for signal events, but can be larger in cc events, (vi) the distance along the beams between 
the reconstructed vertex of the B candidate and the vertex of the other tracks in the event. This 
is consistent with zero for continuum events, but is sensitive to the B lifetime for the signal events. 

The NN is trained with simulated continuum and signal events. We find agreement between 
the distributions of all six variables in simulation and in control samples of off-resonance data and 
of B —> Dtt events. We apply a loose pre-selection on the NN (0.4 < NN < 1.0) with a 90% 
efficiency on signal and a 68% rejection power over continuum, but then use the NN itself in the 
likelihood fit to fully exploit its discriminant power. 

A B candidate is characterized by the energy-substituted mass tties = ^/(f + Po • Pb) 2 /E 2 — p 2 B 
and energy difference AE = E* B — \yfs, where E and p are energy and momentum, the asterisk 
denotes the CM frame, the subscripts and B refer to the initial e + e~ state and B candidate, 
respectively, and s is the square of the CM energy. For signal events mES is centered around the 
B mass with a resolution of about 2.5 MeV/c 2 , and AE is centered at zero with an RMS of 0.017 
GeV. 

Considering both the right sign and the wrong sign sample, 28621 events survive the selection 
described above and the loose requirements \AE\ < 100 MeV and m-ES > 5.2 GeV/c 2 . The dominant 
background still comes from continuum events, but we also need to take into account background 
from Y(4S) — * BB ("BB") events. We consider separately the B — > Dtt background since it 
differs from the signal only in the AE distribution. This decay mode has a very low value of tb 
(~ rB(D°K)\ 2 , where A ~ 0.22 is the sine of the Cabibbo angle), and therefore in the likelihood 
fit we will consider it as a background only for the right sign sample. 



3 Likelihood Fit and Results 

The signal and background yields are extracted by maximizing the extended likelihood C = 
e~ N ' UiLi Ci{xi)/N'\. Here N' = N DK + N cont + N BB + N Dn is the sum of the yields of the 
signal and the three background contributions, af = {AW, AE, m^s}, and the likelihood of the 
individual events (£j) is defined as 



NjyK N^^ 

1 . D /m(^) + i I C °n fconti^) (3) 
1 + rtADS 1 + ticont 



10 



for right sign events and 




NcontRcont ?WS( 
-i _j_ r> Jcont \ 



for wrong sign events. 



We have defined R parameters for the backgrounds with the same definition as in Eq. ^ The 
individual probability density functions (/) are derived from the MC and the three variables are 
considered uncorrelated in all cases, apart from m-ES and AE for the Dir background, since the 
correlations are not negligible. For the latter we have therefore utilized a two dimensional non- 
parametric distribution . The NN distributions are all modeled with a histogram with eight bins 
between 0.4 and 1. The w-es distributions are modeled with a Gaussian in the case of the signal, 
a threshold function JH] in the case of the continuum background, and the sum of a threshold 
function and a Gaussian function with an exponential tail in the case of the BB background. 
Finally, the AE distributions are parametrized with the sum of two Gaussians in the case of the 
signal, an exponential in the case of the continuum background, and a double exponential in the 
case of the BB background. 

We perform the fit by floating the four total yields (Njjk, N cont , Nbb, and Nd w ), the three 
R variables and the shape parameters of the threshold function used to parametrize rn.Es for the 
same and opposite sign continuum background. Figure 121 shows the distributions of the three 
variables in the selected sample (separately for same sign and opposite sign events), with the 
likelihood projections overlaid. The fit yields Rads = 0.012±g:gjg, N DK = (14.7 ± 0.6) x 10 2 , 
N cont = (239.3 ±2.1) x W 2 ,N BB = (25.5 ±1.6) x 10 2 , N Dn = (6.7 ±0.4) x 10 2 , R cont = 3.05 ±0.07, 
Rbb = 0.42 ±0.07. 

Equation |21 assumes that the efficiencies are the same between the right and the wrong sign signal 
samples, regardless of the different Dalitz structure. This has been tested on MC and proved to be 
true within a statistical error of 4%. We then consider this as a systematic error on Rads- We also 
repeated the fit by varying the probability density function parameters obtained with MC within 
their statistical errors and by estimating f^J^ S on off-resonance data and f^ WS on exclusively 
reconstructed Dir events. To account for the observed variations, we assign a 0.008 systematic error 
on Rads- The uncertainty due to B decays with distributions similar to the signal, in particular 
B -» dWtt, D*K, D^K*, and KKtttt , is estimated by varying their branching fractions within 
their known errors and found to be 0.00006 on Rads, and therefore negligible. The absence of 
further modes which might fake signal has been checked comparing data and MC samples in the 
sidebands of the AE and m D o — m^o distributions. 

Following a Bayesian approach, we extract r# by defining the a posteriori probability 



where £ = C00S7, Rads{tBi r D, C) is given in Eq. ^ and p(re,rD,£) is the a priori probability 
for these three quantities. They are considered uncorrelated, with £ and r# distributed flat in the 
range [—1,1] and [0,1] respectively. The a priori probability for rp is a Gaussian consistent with 
rjj = (0.214 ±0.008 ±0.008)% 0. The likelihood C(Rads) is obtained by convolving the likelihood 
returned by the fit with a Gaussian of width 0.0076, equivalent to the systematic uncertainty. 

Figure shows C(Rads) and £(re). We set a 95% confidence level (C.L.) limit, by integrating 
the likelihood starting from Rads = (re = 0), thus excluding unphysical values, and we define 



£(r B ) 



J p{rB,r D ,£)£(R A Ds(rB,rD, Q)dr D d£ 



(4) 



/ p{r D o K +,r D ,t,)C{R A Ds(rB,r D , £))dr D d£dr B 
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the 68% C.L. region, for each variable r = Rads or r B-> as the interval where C{r) > C m i n and 
^% = J c(r)>Cmin C(r)dr. 

4 Conclusions 

In summary, we measure the ratio of the rate for the — > [K^tt^tt^dK^ decay to the favored 
decay B ± -» {K ± ^i: a } D K ± to be Rads = 0.012±^(stat)±[];^(sys). This result is consistent 
in central value and similar in sensitivity with our completely independent previously published 
result fB]. The measurement is not significant and therefore we set a 95% C.L. limit Rads < 0.039. 
We use this information to infer the ratio of the mag nitudes of the B~ D°K~ and B~ -► D°R- 
amplitudes to be re = 0.091 ± 0.059 and consequently set a limit rs < 0.185 at 95% C.L. 
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Figure 2: Likelihood fit projection of the NN, AE, and toes distributions separately for the same (top) and opposite (bottom) sign 
sample. To visually enhance the signal, the distributions for the latter sample are shown after cuts, with a 67% signal efficiency, on 
the ratios between the signal and the signal plus background likelihood of all the variables other than the one shown. The points 
with error bars represent the data while the dashed, dash-dotted, and solid lines represent the contribution from continuum, BB, 
and Dtt background, respectively. The dotted line represents the signal contribution, visible only in the same sign sample. 




Figure 3: Likelihood as a function of Rads (left) and of tb (right). While the left plot shows the 
actual experimental result of the measurements, the right plot is obtained in a bayesian approach 
assuming flat prior distributions for re, C and 7. The 68% and 95% region are shown in dark and 
light shading respectively. 



15 



